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Abstract 

In this paper we present the full theoretical model of a modified Ives-Stillwell experiment where 
counter propagating lasers are used to form a narrow interference fringe when the lasers form a 
double resnonance. This narrow resoannce can be as small as 1 Hz wide and therefore provides a 
connection between the atomic resonance in its rest frame and the laser frequency in the lab frame. 
The current paper builds on a simplified approach suggested recently [l| and presents a fully 
developed theory of the interaction within the Lorentz vilating electrodynamics of the Standard 
Model Extension [2]. 

PACS numbers: 
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INTRODUCTION 



Precision tests of special relativity are becomming more important as we continue to 
discover more about cosmology and particle physics. Dark matter, dark energy for example 
has a growing experimental support and little strong theoretical support. Theories that 
attempt to reconcile them with quantum electrodynamics (QED) and the standard model 
tend to introduce Lorentz violating features. Thus there is a growing suspician that a 
sufficently sensitive experiment will detect a violation. In an attempt to reconcile QED with 
potential violations of Lorentz invariance, a parameterisation of potential Lorentz violating 
terms has been attempted [2| in a Standard Model Extension (SME). The benefit of this 
parameterisation over previous versions is that it allows an experimental analysis to be 
built from Maxwell's equations. Many of the parameters in this formalism can be highly 
constrained by existing experiments, such as astronomical observations and accelerator tests. 
This leaves a few parameters that are both accessible to desktop optical experiments and 
cannot be easily constrained using other methods. The experiments that have access to 
these parameters are the Michelson interferometer, the Kennedy-Thorndike test and the 
Ives-Stillwell test. Of these, the Ives-Stillwell test appears to be the only experimental test 
that is capable of accessing a parameter called Kt r , and for this reason this is the least well 
known parameter of the SME. 

The Ives-Stillwell experiment is a precision measurement of the Doppler shift of light. 
The modern version of this experiment uses counter propagating lasers that are arranged to 
be co-linear with an accelerated ion beam. Tuning the lasers into resonance with the beam 
allows the lab frame frequency to be compared with the atomic rest frame frequency. Modern 
laser Ives-Stilwell experiments using storage-ring facilities have provided the tightest bounds 
on K tr so far with an upper bound of 3 x lO -8 ^. Current predictions for future storage- ring 
experiments suggest that an upper bound of K tr < 1 x 10~ 9 might be achievable in these 
techniques (J). The current paper explores the theory of a new version of the Ives-Stillwell 
experiment using an EIT resonance in a sample of moving atoms as a potential tool for 
further improving this measurement. 

In the framework of the SME, the sensitivity to an Ives-Stilwell experiment is given by 
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the following formula [5j 



veVw 



1 + 2K tr (j3 2 + 2/3 ffi • /?) , 



(1) 



where ve and ^ are the Doppler shifted frequencies of the two counter propagating laser 
fields, equivalent to u a and v p only now defined to propagate along the axes of East and 
West in the laboratory frame: the orientation of the experiment in the laboratory frame is 
significant due to a sensitivity to the Earth's velocity (labelling the fields in this manner 
removes any ambiguity between the direction of propagation and the coupling and probing 
fields required by EIT which will be encountered later). The rest frame transition frequency 
is z/ , n tr is the SME parameter under investigation, v = /3c is the observer's velocity in the 
laboratory frame and /3q is the Earth's velocity in a frame centered upon the sun. It can be 
seen from equation [T] that there are two access channels to n tr for experiments of this type, 
one relying solely on the observer's laboratory velocity j3, the other on the vector product 
of the observer's velocity in the laboratory frame with the laboratories velocity in the Sun's 
frame /3 ffi . The sensitivity of the Ives-Stilwell experiments can be expressed as, 



In these experiments the frequencies Ve, vw and Pq are measured independently, this means 
that any uncertainty in frequency determination enters in quadrature. Although this can 
introduce significant systematic uncertainties, referencing to a frequency comb removes this 
as a limiting factor in these experiments. From equation (jSJ), it is clear that the sensitivity 
of experiments scales quadratically with j3, implying that increased velocities will lead to 
increased resolution. However, there is an element of diminishing returns associated with 
using extremely high observer velocities as relativity itself restricts the velocities attainable 
by virtue of E 2 = p 2 c c + m 2 c 4 . 

At present the sensitivity of Ives-Stilwell experiments is limited by the linewidth of the 
Lamb dip which is 10.8MHz, about three times the natural linewidth of the transition. As 
a result of the relatively large linewidth current experiments are approaching the maximum 
sensitivity attainable by this detection scheme. 

One way to overcome the limitations imposed by relativity on the observers velocity, 
and saturation spectroscopy on the observed linewidth, is to move away from fast beams 
in favour of slower beams that enable one to create narrower linewidths. In this paper 



veVw 



1 + 2n tr p 2 . 



(2) 
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we present a detailed experimental analysis of a new type of Ives-Stilwell apparatus which 
we will call the Modified Ives-Stilwell (MI) experiment. This takes a different approach 
to other experiments which are sensitive to a variation in the speed of light. However the 
concept is similar to the traditional Ives-Stillwell experiment in that a moving observer is 
used to measure the Doppler shifted frequency of two counter propagating laser fields. In 
this experiment a coherent interaction between the atom and the two laser fields creates an 
extremely narrow feature in the absorption spectrum. This can be orders of magnitude less 
than the excited state linewidth, and can be exploited as a precision frequency discriminator 
in a search for lorentz violating Doppler shift. 

For the MI experiment, the optical fields counter propagate with respect to one another 
and the moving observer is a beam of 85 Rb atoms. The two laser fields are chosen to have 
precisely the same laboratory frequency Vq. 
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(a)Laboratory frame: (b)Atomic frame: 

FIG. 1: Experimental premise. In the rest frame of the laboratory the lasers used are far from 
resonance from either transition, however for atoms with an appropriate velocity both lasers are 
Doppler shifted into resonance forming a A-transition. A coherent interaction between the two 
optical fields and moving observer results in a coherent dark resonance in the absorption spectrum 
of one of the optical field which is extremely narrow. The line centre of the resultant coherent 
resonance is sensitive to small frame dependent changes in the speed of light. 



The co-linear alignment of the lasers is maintained by a folded interferometer and the 
laboratory frequency vq is chosen to be halfway between either leg of a A-system (see Figure 
[I]). The A scheme is defined by the F = 2 «-> F' = 2 and F = 3 «-> F' = 2 hyperfine 



levels in the 5 2 Si/2 5 2 Pi/2, D\ line, in 85 Rb. The lab frame laser frequency is Uq 
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377.107, 408, OTH z. The critical velocity is defined by the ground state hyperfine splitting 
for the Di line and is v = 1207ms -1 for 85 Rb and atoms with the critical velocity are trapped 
in a dark state by the combined interaction of the two laser fields. 

The laser can be stabilized using an optical frequency comb and the linewidth of the 
coherent dark resonance we have achieved is less than 300Hz in the current apparatus. Thus 
we are able to compare the lab frame fequency with the atomic frame frequency with high 
accuracy. 



THE PHOTON SECTOR OF THE MINIMAL SME 

In this section we breifly review the SME so that the terminology of the derivation that 
follows is placed in context. Kosteleckyand collaborators construct the Lorentz violating 
extension to QED by considering only CPT and Lorentz violating terms arising from the 
mrelv photon sector of the full SME. From this one arrives at the following Lagrangian 
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2 = - 4 V*" + 2 ^ ^ A F " U ~ 4 (** W F F " r ( 3 ) 
In equation the electromagnetic field tensor is defined conventionally as F^ v = (d^A u — 
dvAfj), while the coefficients (/caf) k and (k F ) K x^ u describe violations of the Lorentz and CPT 
symmetries respectively. In the limit that these coefficients are zero, equation (EJ) reduces 
to the familiar QED Lagrangian in terms of alone. The coefficient Jj&p)^^ describes 
Lorentz violation. It shares all the symmetries of the Riemann tensor [9J and contains 19 



independent parameters [8]. The coefficient (/caf)* is odd under CPT conjugation and has 
dimensions of mass. 

There is an intimate relationship between the Lorentz and CPT symmetries; the con- 
struction of the CPT symmetry is Lorentz invariant. Hence Lorentz violation permits, but 
does not require CPT violation [lo|. Strong constraints have been placed on (/caf) k via 
spectropolarimitary of distant galaxies 0, 1 111 ] and it has also been shown that finite values 
of (Avlf) k make the Lagrangian of equation ([3]) unstable. It is therefore reasonable to ignore 
the CPT violating contributions to equation Consequently, setting (/caf) k = 0, greatly 
simplifies the Lagrangian of equation ([3]), leavin a purely Lorentz violating contribution: 

4 ^ 4 



& = - 7 V" - 7 (** W FKXP " V - ^ 
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A physical interpretation of the tensor (k E ) KXfll/ in equation (j4]) containing non-zero terms is 
not obvious. However, there is a particularly useful decomposition of (k E ) KXflu which can be 
used to present an analogy between Kostelecky's Lorentz violating extension to QED and 
the conventional condition of light propagating in homogeneous anisotropic media 



KDE 


,jk _ 


-2(k E ) 0j0k , 


K-HB 


J k = 




KDB, 


jk _ 


(k F )°^e k pq , 




,jk _ 





(5) 

From these 3x3 matrices (kde)^, (i^hb)^, {^db)^ and {khe)^ a relationship between the 
electric displacement and auxiliary magnetic fields D and H, and the electric and magnetic 
fields E and B can be derived via elements of the Lorentz violating parameter {k E y x ^ 1 ' 



1 + I^DE KDB 
KHE 1 + Khb 



)l 


: 







(6) 



The Lagrangian of (jlj) can then be expressed in terms of only the Lorentz violating param- 
eters of equation (jSJ) and the electric and magnetic field components as §}, 

& = l - (E 2 - B 2 ) + ^E • (k de ) • E - • (k HB ) B + E (k db ) • B. (7) 

A further decomposition of {k E ) K x^ v defines a set of parameters which are directly accessible 
via experimentation {§]], 

{K e+ ) ik = ^DE + K HB ) jk , 

(n e _y k = 1 -{k de - n HB y k - 1 -&\k de ) 1 \ 

(K 0+ ) jk = -(K DB + K HE ) ik , 
{K -) ik = ^{.KDB - K HE y k , 

«tr = ^De) U - (8) 

Each of the K-parameters defined in equation (jBJ) represents a 3 x 3 matrix except n tr which 
is a scalar quantity. Parity odd terms are contained within the antisymmetric k q+ and 
symmetric k _ and the parity even terms in the symmetric K e+ , n e _ and K tr . Using the 



definitions given in equations ([6]) and equations dHJ) it can be shown that the Lagrangian in 
equation (j3j) can be re-expressed in such a way so as to make a direct connection between 
Lorentz violating parameters and experiments Q], 



Inspection of equation (jHJ) immediately indicates a relationship between n tr and the effective 
permittivity e and permeability fi, 



Equation ffTUj) can, albeit via an enormous oversimplification of the problem, be related to 
the one-way speed of light, 



Equation [TT] illustrates how K tr acts somewhat like a refractive index, though in order to 
describe u correctly requires an appropriate definition of reference frame in which u is defined. 

The one-way speed of light 

The experiment described here, as well as conventional Ives-Stilwell experiments, are lead- 
ing order sensitive to deviations in the one-way speed of light from the spacetime constant 
c. Although equation (1111) illustrates a connection between K tr and the speed of light, with 
K tr acting somewhat like a refractive index. In order to corectly describe the experiment a 
full derivation in a Lorentz violating electrodynamics is required. This has been performed 
in a detailed analysis starting from the SME. 

Using the definitions of the Lorentz violating electric displacement vector and auxiliary 
magnetic fields from equation (jH]) one can arrive at a Maxwell's equations, using the modified 
D and H fields from equation [6], 




(9) 




(10) 



d s Bi 

e^dM - d Di 
e ikl d k B x + doB* 



0. 



o. 



(12) 
(13) 
(14) 
(15) 
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From these, a modified dispersion relation for the electric field can be derived using the 
Lorentz violating analogue of Amperes law, equation 



plane-wave with four-momentum p a = (p°,pi), [a. Il2l|: 



14]) . Therefore, for an electromagnetic 



M jk E k = (_ 5 iy _ pjp k _ 2(k F y^ k Pl3Pl )E k = 0. (16) 

A condition of equation (TIB"]) is that the determinant of the coefficient should vanish, |M J | = 
0. Solving \M jk \ = one arrives at a relation between p° and pP |8[], 

p° ± = (l + p±a)\^l (17) 

where p = — ^k a a , a 2 = \{k a p) 2 — p 2 and with, 

k * = (k F ) a ^p,p u , (18) 

and p^ = p^/lpl- From this dispersion relation the phase velocity of light, defined as u = 
PoP* Iv 1 1 can then be determined, 

u = c(l + p±a). (19) 

This expression is important. It shows how, in this particular framework, deviations of the 
phase speed of light, u, from the spacetime constant c are contained within the parameters 
p and cr. The ability to describe an optical experiment in terms of the SME is dependent 
on evaluating both p and a in terms of parameters described in equation ([8]); this requires 
a determination of k a @ in terms of these parameters also. Although k a/3 must have been 
calculated by others previously, an explicit expression for it has never been published and 
so one is presented here. 

The tensor k a/3 in equation ( Tl8l) is a function of the wave unit-four-momentum, a conse- 
quence of which is that the terms p and a from equation ffl9|) are both frame dependent. 
Therefore a meaningful evaluation of either p or a, and subsequently u, requires the specifi- 
cation of a reference frame. Although this reference frame can be freely chosen, calculations 
can be significantly simplified by choosing a special frame in which to evaluate k a/3 and 
then making the appropriate transformation into an agreed frame where experiments are 
analysed and their results compared. For the case of an Ives-Stilwell experiment a suitable 
frame to choose is one at rest in the laboratory frame with it's third spatially axis aligned 
co-linearly with the two optical fields. In this special frame, which shall be refered to as 
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East West 
p'» = (1; 0, 0, - *) ) p'" = (1; 0, 0, 1) 



Vat 



FIG. 2: The Experimental frame, Optical fields propagate with phase speeds ue and u\y along 
the unit four-vectors p !>M = (1; 0, 0, —1) and ]5' M = (1; 0, 0, 1) respectively. The Easterly and Westerly 
unit vectors in a the laboratory frame line up with p'^ = (1; 0, 0, —1) and p'^ = (1; 0, 0, 1). 

the experimental frame and denoted by <f , the wave unit-four-vectors of the counter prop- 
agating fields therefore take the form p ftl = (1; 0, 0, ±1), see Figure [2J Note that the third 
spatial component of p >fl to be +1 for fields propagating from East to West, and —1 for fields 
propagating from West to East throughout. 

The initial aim of this derivation was to be able to describe the phase speeds ue and uw, 
for the two counter-propagating lasers used in this experiment in terms of the parameters 
K e+ , K e _, n 0+ , k _ and K tr , defined in equation ([S]), as this would then allow easy comparison 
of the results from this experiment with complimentary experiments. 

However, it is considerably easier and generally more useful to apply frame transforma- 
tions to the 'electrodynamic' parameters kde, k-hb, ^db and khe defined in equation ([5]). 
These parameters are explicitly defined in terms of (kF) afJ-f3u and are therefore more easily 
related to both p and a, a conversion into K e+ , K e _, k q+ , k q _ and K tr via equation ([8]) then 
becomes possible. 

We begin the determination of p and a by rearranging equations §5§ to give (kp) ^^ 
terms of the matrices kde, &hb, ^db and khe'- 

(k F ) ojok = -\{K DE y\ 



(k F r rs = -e k rs e/ q (K HB y\ 

{k F f m = \{K DB y k e k pq . (20) 

It should be noted that in order to completely describe hp using these equations it is 
extremely important to remember it has the same symmetry properties as the Riemann 



9 



tensor 



(k F ) a ^ = -{k F f aixv = -(k F ) a ^ = (k F )^ (21) 

and 



{k F ) v + {k F ) a ^ + (k F ) av ^ = (22) 

Folowing this it is possible to fill all non-zero terms of in (k F ) a ^ tiU with elements of 
kde, khb, i^db and Khe- The tensor (k F ) a/3 ^ u could then be contracted with an appropriate 
wave unit-four- vector according to equation (|18|) to leave k al3 for both Easterly and Westerly 
propagating fields also in terms of kde, ^hb, ^db and khe 3]- From which the Lorentz 
violating contributions to u F and uw can be determined. 

With an expression for k a ^ available in terms of kd F , khb, k db and Khe alone, and after 
evaluating p in terms of k af3 , 

p = -\K a = 

= - l - [k 00 - k 11 - k 22 - k 33 ] , (23) 

it was possible to rewrite equation (1231) as for a beam propagating along the Easterly unit 
vector, p^ = (1; 0, 0, 1), as: 

P=-\ [("de - khb) U + {kde - khb) 22 + 2{k db ) 21 - 2{k db ) 12 ] . (24) 

To evaluate p for p^ = (1; 0, 0, —1) simply change the sign of the k^B terms in equation ( 1241) . 
Similarly, by evaluating a in terms of k a @, 

o 2 = \{K,) 2 -p 2 

= [(k 00 ) 2 + (A; 11 ) 2 + (k 22 ) 2 + (k 33 ) 2 

+ (jfcl2)2 + (A ,13 )2 + (fc 21 )2 + (fc23) 2 + (fc 31)2 + ^32)2 
_ ( fc 01)2 _ (A; 2)2 _ (A ,03 )2 _ (A ,10 )2 _ (A;2 )2 _ (A ,30 )2] 

- l -[k™-k ll -k 22 -k 33 ] 2 , (25) 
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and evaluating in terms of kde, ^hb, ^db and khe for = (1; 0, 0, 1), 

a 2 = 1(4[(A; D b) 11 ] 2 + 4[(A; DB ) 12 ] 2 + 8(k DB ) 12 (k DB ) 21 + A[{k DB ) 21 } 2 - 8(k DB ) u (k DB ) 22 
lb 

+ 4[(^b) 22 ] 2 - ±(k DB ) 12 (k DE ) u - 4(k DB ) 21 (k DE ) u + [{k DE ) n f + S^b) 11 ^) 21 

- 8{k DB ) 22 {k DE ) 21 + ±[{k DE ) 21 ] 2 + ±{k DB ) 12 {k DE ) 22 + A{k DB ) 2 \k DE ) 22 

- 2{k DE Y\k DE ) 22 + \{k DE ) 22 ] 2 - 4(k DB ) 12 (k HB ) n - 4(k DB ) 21 (k HB ) n 
+ 2(k DE ) u (k HB ) n - 2(k DE ) 22 (k HB ) n + [{k HB ) n \ 2 + 8(A; Di? ) 11 (A;^) 21 

- 8(k DB ) 22 (k HB ) 21 + 8(k DE ) 21 (k HB ) 21 + A\{k HB ) 21 } 2 + 4(k DB ) 12 (k HB ) 22 

+ A(k DB ) 21 (k HB ) 22 - 2(k DE ) u (k HB ) 22 + 2(k DE ) 22 (k HB ) 22 - 2(k HB ) 11 (k HB ) 22 + [(k HB ) 22 ] 2 ). 

(26) 

Evaluating a is considerably more effort than p, but thankfully, it can be neglected for the 
purposes of this experiment. 

The tensors k a/3 can be compared to the dielectric tensor of an anisotropic crystal from 
classical electrodynamics. Equations (1231) and (|25l) show, p and a are related to the trace 
and off-diagonal terms of k al3 , which can be related to the phase speed and isotropy, and 
birefringent contributions from Lorentz violation respectively. 

EIT IN A BEAM OF 85 RB ATOMS 

The MI experiment seeks to measure any small deviations in the speed of light which 
manifest themselves as an asymmetric Doppler shift between two counter-propagating laser 
fields. In order to make any statement about the validity of the Doppler shift predicted by 
special relativity, knowledge of the observers precise velocity is required. 

A highly idealized system 

Initially consider a beam of idealized three level atoms travelling with a velocity v c and 
interacting with two counter-propagating laser fields of equal rest frame frequency vq and 
orthogonal circular polarisations \cr±). The frequency z/ is chosen such that the approaching 
field will be blue shifted into resonance with the z/ 31 transition while the receeding field will 
be red shifted into resonance with the z/32 transition. For a symmetric Doppler shift of the 
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m F (3),g F (3)) 




m F (l),5 F (l)) 



(b)The laser fields form a A-system 
with the three level atom 



(a)Two counter-propagating lasers 
interact with a three level atom 
travelling co-linearly with velocity v c 

FIG. 3: EIT in idealized three level beam. The lasers are assumed to have equal rest frame 
frequency uq chosen such that for an atomic velocity v c the weak probe field and strong coupling 
field with orthogonal circular polarisations form a A-system with the three energy levels of the 
atom. 



two fields this occurs for a laser frequency of, 

^31 + ^32 

in an atom travelling with velocity v c defined as, 



(27) 



^21 



^31 — ^32 



-C. 



(28) 



2z/ u 31 + z/ 32 

If the states |1), |2) and |3) of our moving observer are composed of Zeeman substates 
with magnetic quantum numbers m F (j) and gyromagnetic ratios g F (j) for state \j) and 
\m F (3) — m F (2)\ = \m F (3) — m F (l)\ = 1 and \m F (2) — m F (l)\ = 2 then the three levels 
can form a A system using orthogonal circular polarisations a + and cr_ , see Figure [31 In 
this model the field co-propagating with the atomic beam is strong and therefore considered 
a coupling field, while the field counter-propagating with the atomic beam is considered a 
probing field. Zeeman sub levels are sensitive to magnetic fields. The change in frequency, 



Au B , of such a state induced by an applied magnetic field is described by 

Au b = ^-m F g F B, 
n 



(29) 
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where B is the magnetic field along the axis of the interaction. An externally applied 
magnetic field can therefore contribute to the detunings Ai and Delta-i- 

For an exquisitely defined atomic beam, i.e. an atomic beam with a velocity distribution 
approximating a delta function, 5(v — v c ), and considering only the first order contribution 
to the Doppler shift, the detunings A and 5 are therefore described by the formulae: 



A(B, v c ) = ia,i (0) + ^ [m F {3)g F {3) - m F {l)g F {\)} B - u {l + v c /c), 
5(B, v c ) = i/ 21 (0) + ^ [mj,(2)^(2) - m F (l)^ F (l)] 5 - 2z/ v c /c, 



(30) 



where ^fe(O) represents the frequency splitting between \j) <-> for 5 = 0. These de- 
tunings provides a description of the absorption as a function of applied magnetic field. 
Figure H] shows the effects of Zeeman detunings in an exquisitely defined atomic beam as 



o 





O 

22 




FIG. 4: The normalized EIT imaginary susceptibility for a hypothetical atomic beam. The beam 
is exquisitely defined with velocity v c and the coupling and probing fields are counter-propagating 
with respect to one another. See text for a description of the parameters used in this simulation. 



described previously. The parameters used are based on a three level system with the 
same hyperfine structure as 85 Rb: u 31 = 377,108,911.7MHz, u 32 = 377,105,876.0MHz 

151 ]. The states are |1) = \F 



and z/ 2 i = 3035.732MHz [15[. The states are |1) = \F = 2 } g F = -1/9, m F = 2), 
\2) = \F = 3,g F = 1/9, mj? = 0) and |3) = \F' = 2,g F = -1/3, m F = 1). The de- 
phasing terms used were 731 = 1 x 10 11 Hz and 721 = 0Hz and the coupling field took a value 
of O.OOOI731. The velocity v c = (z/ 31 — ^32)/ (^31 + ^32) ~ 1207ms -1 . 
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Doppler broadening 



An exquisitely denned velocity distribution is of course highly unrealistic as the velocity 
distribution will governed by Maxwell-Boltzmann statistics. Subsequently, there will be a 
large number of atoms of varying velocities for the coupling and probing fields to interact 
with. With the Doppler width of rubidium being some hundreds of MHz one might expect 
the narrow transmission window depicted in Figure H] to be washed out. However, numerical 
simulations have shown this not to be the case. 

The Maxwell-Boltzmann velocity vector distribution function is described by, 

This describes the probability, P(vj), of an atom travelling along the unit vector j with 
a speed Vj, T is the absolute temperature, k is the Boltzmann constant and m is the 
atomic mass. In the cell there will be a number of atoms of different velocities and therefore 
detunings, A and S, contributing to the overall observed susceptibility. This can be simulated 
by integrating the product of the imaginary part of the susceptibility with the velocity 
distribution function, 

/ = jQfcMfrvjjlPMdVj. (32) 

The critical velocity v c = 1207ms -1 for the D\ line of rubidium is far out in the tail of the 
Gaussian velocity distribution where the distribution is quite flat. Because both lasers are 
derived from the same source, and therefore have the same Uo, an atom with velocity v^v c 
will see an equal and opposite Doppler shift away from resonance for the two optical fields. 
Coherent phenomena such as EIT or coherent population trapping have previously been 



shown to exhibit velocity selective behaviour 161 ] . A simple argument to explain the velocity 



selective nature of EIT in our scheme is depicted in Figure El For an atomic beam travelling 
with velocity v c the two photon Raman resonance condition is satisfied and EIT can occur. 
However, for small deviations from this critical velocity the Raman resonance is not fulfilled 
by virtue of opposite Doppler shifts in the two interacting beams, subsequently EIT cannot 
be generated. 

In addition to this qualitative argument, a numerical integration of equation (1321) was 
performed for one hundred points ranging from —40 — > 40nT between the limits of ±4.4ms _1 
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, see Figure [H These limits are denned by the excited state linewidth which is 5.746MHz. 
The lineshape which results from this numerical integration is indistinguishable from that of 
an exquisitely defined atomic beam of velocity v c . It is as though the velocity distribution 
has collapsed into a delta function making the transition velocity selective; only those atoms 
with a velocity close to v c contribute to the observed signal. For a bariety of conditions, the 
resultant lineshape remained unaffected for integration limits up to several hundred MHz. 
The transition is therefore extremely insensitive to Doppler broadening. 

THE IVES-STILWELL EXPERIMENT AND THE SME 

To analyse the MI experiment in the SME we begin by defining the Doppler shifted 
frequencies of the two counter-propagating fields in the reference frame of the atom. For an 
experiment as depicted in Figure [2] and including any modifications to the phase velocity 
which might occur from a violation of Lorentz invariance, i.e. ue,w, the Doppler shift can 
be described as, Where the phase velocity of light is given by equation [19] and v F and vw 
refers to the frequency of light seen by an observer traveling along an Easterly and Westerly 
unit vector respectively, with an observer also propagating along the westerly unit vector. 

The Modified Ives-Stilwell experiment 

In equations [30] we showed that the one and two photon detunings, A and 5, can be 
described by, 

A(B, v c ) = u 31 (0) + ^ [m F (3)g F (3) - m F {l)g F {l)\ B - u (l + v c /c), 
6{B,v c ) = z/ 21 (0) + ^ [m F {2)g F {2) - m F {l)g F {l)\ B - 2z/ v c /c. (33) 

However, these detunings do not take into account any Lorentz violating contributions. In 
order to do so, the spacetime constant c appearing in equation (!33l . which arrises from 
the Doppler shift, must be replaced by the appropriate phase speed u F ,w from equation 
( I19p . For an atomic beam propagating Easterly with a probe field propagating Westerly, 
and expanding only to leading order in atomic velocity f3 at and p, the one and two photon 
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detunings A and 5 can be written as: 

A(B,p at ) = z/ 31 (0) + ^ [mjr(3)^(3) - m F (l)g F (l)] B - u {l + f3 at ) + VoPatPE, 

= U 31 (B) - Uq(1 + (3 at ) + VOPatpE, (34) 



6(B,p at ) = z/ 21 (0) + ^ [m F {2)g F {2) - m F {l)g F {l)] B - 2v (5 at + v (3 at {p E + p w ), 
= v 2X {B) - 2v f3 at + u p at (p E + p w ), 

= v 2 i(B) - 2v (3 at - 2K tr v (3 a t, (35) 

where f3 a t = v at /c. The final terms in equation ( 13~4"1) and ( 1351) contain the Lorentz violating 
contributions to the detunings. As one would hope, in the limit that K tr = these detunings 
reduce to those predicted by special relativity. 

Simulations of the expected EIT signal for finite values of n tr revealed it to be far less 
sensitive to a Lorentz violating contribution to the single photon detunings, A, than it is 
two the photon detunings, 5. These simulations demonstrated that the effect of even large 
values of K tr in this context had little effect on the observed spectra, with noticeable effects 
for the current experimental resolution only arising for values of Kt r orders of magnitude 



larger than the current upper bound, K tr < 8.3 x 10~ 8 [171 ]. 

In contrast, the contribution of K tr to the observed EIT spectra from the two-photon 
detuning 5 was very different, having an enormous effect on the magnetic field at which the 
peak of the spectrum occured. This effect persisted even after integration over the thermal 
velocity distribution, confirming that the resonance remains velocity selective even in the 
presence of Lorentz violation. 

From equation fl35l) it can be seen that the offset of the peak, 5, arrising solely from 
Lorentz violating effects, for such an experiment is: 

C = 2K tr Vol3at, 

~3x l^K tr Hz. (36) 



For the most magnetically sensitive transition in the Dl line of 85 Rb an external magnetic 
field will cause a change in the two photon detuning of ~ lMHz/Gs = 10 w Hz/Tesla, so 
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the offset described by equation (136|) is approximately equivalent to: 



2Kt r Vof3at, 



0.3K tr Tesla. 



(37) 



For a value of n tr consistent with the current experimental upper bound this would cor- 
respond to a shift of 25nT in the peak of an EIT signal. Magnetometry at this level of 
sensitivity is already achievable, with magnetometers relying on phenomona such as EIT 
being shown to exhibit sensitivities already on the order of 10pT/ \/ Hz frsl. \\\ 



Although an obvious experiment to measure K tr might be to look for the deviations of an 
EIT resonance line centre from OT magnetic field in an atomic beam. There are a number 
of factors, other than Lorentz violation, that could cause a shift in the line centre of such 
a transition, with stray magnetic fields being of particular concern. Such a measurement 
would be fraught with systematic uncertainties. An alternative, which is capable of removing 
many of these systematic effects, is to make a differential measurement. 

For the model considered here, there should be two atomic beams propagating in op- 
posite directions along the experimental axis which could support an EIT resonance. This 
system can be modeled by two identical experiments with the probe and coupling field de- 
tunings interchanged between them, see Figure [HJ To begin the analysis we first consider 
a toy-model for this system which consists of two similar A-atoms, one propagating East 
and one propagating West, each with speed v c . Assume that the A-atoms are identical 
except that the ground states states |1) and |2) switch roles between the true ground, and 
metastable ground states and that the appropriate transitions are addressed by the appro- 
priate polarisations for each laser, see Figure M A finite value of n tr manifests itself as an 
atomic-beam direction-dependent magnetic field, with the Lorentz violating contribution to 
the two-photon detunings given by: 



where the superscripts on 5 represent the propagation direction of the atomic beam. A 
comparison between the line centres of two such systems could therefore be used to perform 
a differential measurement where there separation could be used to measure the size of Kt r 




= -2u j3 at K tr 



(38) 
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according to, 



AH 



AK tr iy (3 at . 



(39) 



Looking for a separation in this way, rather than an absolute offset from zero, removes many 
of the systematic contributions to the uncertainty in the measurement. 

In a real atomic system the direction of the shift, i.e. whether the line center moves to 
a positive of negative magnetic field, is also determined by the internal structure of each 
A-transition contributing to the system. 

TOWARDS A REAL ATOMIC SYSTEM 

So far we have described how EIT can be generated in an idealized atomic beam selected 
from a thermal distribution, however this is far from the situation encountered in a real 
atomic system. Before attempting to tackle the D\ line in 85 Rb (i.e. the system used in 
this experiment) we considere the somewhat simpler case of two A-systems formed between 
magnetic sublevels interacting with the same optical fields Q p and Q c as illustrated in Figure 
[TU1 In this system the energy levels |1) and |4), |2) and |5), |3) and |6) are degenerate with 
one another, except for some small detuning which can arise from an externally applied 
magnetic field. The Hamiltonian of the combined system can therefore be written as, 



where the diagonal terms are 3x3 matrices with the appropriate choice of detunings and 
Rabi fields for the respective A-system inserted. To examine this system in more detail 
we construct a Master equation using all six energy levels making sure to also include an 
appropriately modified decay operator, 



Dp = -[2r 31 p 33 |i)(i| + 2r 32 p 3 3|2)(2| - 2(r 31 + r 32 ) P33 |3)(3| 

+ 2r 3lP66 |4)(4| + 2r 32 p 66 |5>(5| - 2(r 31 + r 32 )p 66 |6)(6| 

- 721 (p 21 |2)(l|+ P 5l|5)(l|+P42|4)(l|) 

-73i(p3i|3)(l| + P6i|6>(l|+P4 2 |4)(3|) 
- 73 2 (p32|3)(2|+p 62 |6)(2|+p4 2 |5)(3|) 

-74iP4i|4)(l| - 752P52|5>(2|- 763 p 63 |6>(3|]+#. C . (41) 




(40) 



1 
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where we have assumed that the two ground states |1) and |4) to remain equally populated, 
pu ~ p 44 ~ 1/2; a reasonable assumption for such a system inside an atomic vapour cell 
where the population is thermally distributed between the two ground states. Therefore the 
general density matrix for this system can be simplified slightly to become, 

/ (l/2 Pl2 Pl3 \ 

P21 p 2 3 



9 



^ Pl4 Pl5 Pl6^ ^ 



V 



\PZ\ P32 J 

Pil P42 P43^ 

P51 P52 P53 

yPei P62 P63 / 



P24 P25 P26 
\^P34 P35 P36 J 

'1/2 p 45 p 46 ^ 

P54 p 56 
\^P64 P65 J 



(42) 



/ 



Solving the Master equation for this density operator we find that the time evolution of the 
coherences occurring in the block wise diagonal terms depend only upon the population of 
their respective quadrant. That is to say that p 31 = /(pu) 7^ f(pu), P64 = /(P44) 7^ /(pu) 
etc. Additionally, because |1) and |4) are connected to |3), and |2) and |5) are connected 
to 1 6) via different spontaneous emission polarisation channels no coherence can be formed 
between them under normal operation. Because the time evolution of the blockwise, off 
diagonal terms, in equation (T42|) does not depend on the level populations, if there is no 
coherence in the first instance, then there can be not be coherence at any time there after. 
This permits further simplification of the density operator, 



/ U/2 Pl2 p 13 ^ 



9 



\ 



P21 p 23 

\^P31 P32 J 







^000^ \ 


\0 Oy 

^1/2 p 45 p 46 ^ 



(43) 



P54 p 56 
yu u uy \^p 64 p 65 J J 

Because the elements of the blockwise diagonal quadrants in equation H3] do not depend on 
any terms contained within the other diagonal quadrant, or the off-diagonal quadrants the 
two systems A and B can be considered independent of one another. Therefore a general 
solution for the steady state density matrix of the combined system is, 



Pss = PaPss + p BPss '» 



(44) 
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where Pa and Pb are the probabilities of being in systems A and B respectively and p ss ' 
and p^g are the steady state density matrices for the two systems if they were isolated. 
Assuming equal probabilities Pa and P B the effective susceptibility of the two systems is, 

2N 

X (1) = B T/ (^13031 + A146P64 + A123P32 + ^&P&b) (45) 

and therefore the effective susceptibility of this combined system is simply that of the sum 
of the two systems A and B. In the D\ line of 85 Rb there are four possible A-systems which 
can be formed using circularly polarised light. We assume that these four systems can 
be treated independently and that the total susceptibility observed will be the sum of the 
susceptibilities of the independent systems. 

3 

EIT in the D 1 line of 85 Rb vapour 

The Maxwell-Boltzmann velocity vector distribution is symmetrical in velocity, therefore 
an atom travelling along the Easterly unit vector in the laboratory frame with a velocity v c 
is as equally likely as an atom travelling along a Westerly unit vector at the same speed, for 
example. Because the rest frame laser frequency has been chosen such that a symmetrical 
Doppler shift brings the two lasers into resonance with their respective transition there will 
be two effective beams which can be selected from the thermal distribution. The population 
inside the cell is initially thermally distributed, this means there are equal populations of 
atoms with their valance electrons in each one of the possible mp magnetic sub levels within 
the F = 2 and F = 3 manifolds of the ground state. The description previously given to 
describe the susceptibility of an idealized three level system assumed that the population 
was initially all in the absolute ground state. Although this is not really the case, only those 
atoms with their valance electron in an appropriate hyperfine state for their orientation with 
respect to the laser fields will contribute to the observed EIT signal. All other atoms can 
only contribute to the background absorption, which changes only on a frequency scale much 
larger than that of EIT. Any optical pumping occurring will only increase the population's 
capable of exhibiting EIT. 
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The separation between the F = 3 and F = 2 states in the 5 2 S"i/2 manifold is ~ 360MHz 
which is large enough to individually probe either of these states without influence from the 
other. Considering only the F = 2 excited state there should be 8 different A-systems which 
can contribute to any observed EIT spectrum; 4 in each of the two atomic beams propagating 
along the experimental axis, see Figure [TT1 Although each A-system interacts with the same 
optical fields, due to different dipole matrix elements as well as different gryomagnetic ratios 
and magnetic quantum numbers each system can support a different linewidth. Regardless 
of the contributing A-system, the line center of each resultant spectrum should be centred 
at zero magnetic field. 

Simulating a violation 

The experiment uses two beams of 85 Rb atoms propagating in opposite directions along 
the East<->West experimental axis. The predicted EIT lineshapes for all eight of the A- 
systems which can be supported by the |1) = \F — 2), |2) = \F — 3) and |3) = \F — 2) 
states within the Dl manifold of 85 Rb have been modelled. Of these eight A-systems only 
six can contribute to any EIT signal, because the internal structure is such that there is no 
magnetic field which can generate a detuning for those systems with |3) = \mp = 0). Of the 
remaining six contributing systems, see Figure [121 only four provide unique deviations of 
line centre from zero magnetic field: For finite values of Ktr transitions with |3) = \mp = 1), 
are equivalent to one another in counter-propagating atomic beams. This also applies to 
transitions with |3) = \mp = — 1). 

The four remaining unique spectra which occur for K tr ^ are shown in Figure [131 The 
linewidths of the induvidual spectra in Figure [T3l differ due to their varying Clebsch-Gordon 
(CG) coefficients, magnetic quantum numbers and gyromagnetic ratios. The A-systems (T) 
and (2) are twice as sensitive to magnetic fields as (3) and @. The simulation assumed a 
value of K tr = 8 x 10~ 8 and a fixed intensity coupling field from which the appropriate CG 
coefficient were derived. All frequency parameters used in this simulation correspond to the 
Dl line of 85 Rb. 

The results of this simulation are presented in Figure [T3l This figure shows each of the 
four unique spectra, their colours are matches to the appropriate A-systems shown in Figure 
[T2l It also shows the summation of all these spectra (black line) taking into account the 
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factor of two contribution in the amplitude of the most magnetically insensitive transitions 
(3) and (4). What this figure shows is that for large enough values of K tr the resultant 
spectrum from the four contributing unique spectra also exhibits a splitting. It also shows 
that for this particular system the separation between resultant peaks in the overall spectrum 
is governed by the splitting between the most magnetically sensitive A-systems (T) and (2). 

Despite the added complexity of this system over the toy-model presented in Figure [9] 
it is still possible to access K tr according to equation fl39l) where the splitting AS is now 
determined by the splitting of the most magnetically sensitive A-systems (T) and (2). 

The dominant contribution to any Lorentz violating signal observed in this experiment 
will originate from the SME parameter K tr . Although there will be additional contributions 
from the other SME parameter, n 0+ , K e _, k _ and n e+ which have so far been neglected 
because of the significantly tighter bounds placed on them by other experiments. 

In this paper we have presented a detailed examination of the Ives-Stilwell experiment 
in the SME for a real atomic system. It shows that the Modified Ives-Stilwell experiment 
can provide a sensitive test of Lorentz Invariance. We have shown that for a relatively 
simple system it is possible to exceed the best current experimental test for n t r. More 
advanced setups would allow considerably better measurements, for example using counter 
propagating atomic beams, rather than a gas cell would enable a further 100 fold increase 
in sensitvity. 
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(a)Exactly on resonance: The atom has precisely the critical 
velocity v c such that the parallel and antiparallel laser fields are 
Doppler shifted into resonance with the transitions 1/32 and ^31 

respectively. 



V32 



> 



< 



1'31 



(b) Small velocity detuning: The atomic velocity is slightly 
detuned from v c such that the resonance frequencies ^32 falls on 
the red detuned side of the laser/excited state linewidth and the 
resonance frequencies P31 falls on the blue detuned side of the 
laser /excited state linewidth. 



> 



V32 



V31 



< 



(c)Large velocity detuning: For atomic velocities far from the 

critical velocity v c the Doppler shift is so great that neither 

frequency is anywhere near the resonance frequencies 1^32 or ^ 31 . 
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FIG. 5: Why this system is velocity selective. Both lasers have the same rest frame frequency vq , so 




-40 -20 20 40 

B/nT 

FIG. 6: Numerical integration of Doppler broadening. This figure compares the imaginary sus- 
ceptibility for that of an exquisitely defined atomic beam (red) and a thermal distribution after 
numerical integration of the atomic velocity from v c ± 4.4ms -1 (black). The integration limit was 
defined by the excited state linewidth of 5.746MHz. Both spectra have been normalized to their 
own maximum value. See text for further details. 
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FIG. 7: A finite value of Ktr shifts the EIT line centre from zero magnetic field to a value S — 
2Kt r vof3 a t- The blue line represents the spectrum in the absence of Lorentz violation, the red line 
represents the spectrum for a finite value of Kt r - 

z 




FIG. 8: The two systems generated in the cell. The experimental axis runs East <-> West. Red and 
blue arrows are the red and blue detuned laser fields respectively. 
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FIG. 9: Toy- model of two counter propagating A-systems. Under the influence of a finite value of 
Kt r the resonance condition of these two systems shifts in opposite directions. 
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FIG. 10: Two A-systems A and B formed between Zeeman sublevels |1), |2) and |3), and |4), |5) 
and 1 6) respectively. The states |1) and |4), |2) and |5), |3) and |6) are degenerate with one another. 
The excited state |3) can spontaneously decay (represented by curly arrows) to either one of the 
two ground states |1) or |4) at a rate of 731 or to the meta-stable state |2) at a rate of I^, while 
the excited state |6) can decay to either one of the meta-stable states at a rate of I"^, or one the 
ground states |4) at a rate of 1^31 
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(a) Thermal population 

F = 3 




(b) Contributing population 

FIG. 11: The population in the cell, is initially thermally populated between the F = 2 and F = 3 
ground states of the 5 2 S±/2 manifold. However, only those atoms with their valance electron in an 
appropriate state will contribute to any observed EIT signal with any other atoms effecting only 
the background absorption which changes only on a frequency scale much larger than that of EIT. 
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FIG. 12: The six A-systems which can be generated in the two atomic beams inside the cell. The 
colour coding of the A-systems corresponds to the spectra depicted in Figure [T3l The colour coding 
corresponds to the predicted spectra that can be found in Figure [13] 
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FIG. 13: Simulating a violation. Displayed are the four contributing spectra. The atomic param- 
eters used in these simulations were matched to those obtained by fitting to current experimental 
data, where we have assumed a Lorentz violation of Kt r = 8 x 10~ 8 (the current best experimental 
limit). The (• • • •) and ( > ) represent transitions between |3) = \mp = ±1) in both velocity 
classes. Because of the symmetrical structure of these transitions there is an additional factor of 
two contribution from these states. The (□□□□) and (+ + + +) represent the |3) = {mp = 2) 
transition in the counter propagating atomic beams. The solid black line is the summed contri- 
bution from all relevant transitions and is the expected output from the experiment assuming an 
equal contribution from all A-configurations. 
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